Quantum Hall Effects in Silicene 
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We investigate quantum Hall effects in silicene by applying electric field E z parallel to magnetic field. Sil- 
icene is a monolayer of silicon atoms forming a two-dimensional honeycomb lattice, and shares almost every 
remarkable property with graphene. A new feature is its buckled structure, due to which the band structure can 
be controlled externally by changing E z . The low energy physics of silicene is described by massive Dirac 
fermions, where the mass is a function of E z and becomes zero at the critical field E a . We show that there are 
no zero energy states due to the Dirac mass term except at the critical electric field E a . Furthermore it is shown 
that the 4-fold degenerate zero-energy states are completely resolved even without considering Coulomb inter- 
actions. These features are highly contrasted with those in graphene, demonstrating that silicene has a richer 
structure. The prominent feature is that, by applying the electric field, we can control the valley degeneracy. As 
a function of E z , Hall plateaux appear at the filling factors v = 0, ±1, ±2, ±3, • ■ • except for the points where 
level crossings occur. 



I. INTRODUCTION 

Electrons in graphene may be viewed as massless Dirac 
fermions in the 2-dimensional space 1 . The study of graphene 
started with revealing an unusual quantum Hall (QH) effect 
under magnetic field 2 9 , where Hall plateaux develop at a se- 
ries of filling factors v — ±2, ±4, ±6, • • • under weak mag- 
netic field, reflecting the existence of the zero-energy states 
and the 4-fold degeneracy associated with the spin symmetry 
and the valley symmetry. Though the 2-fold spin degeneracy 
may be resolved in a magnetization experiment, there is no 
way to access to the valley symmetry in graphene. 

Recently silicene, a monolayer of silicon atoms form- 
ing a two-dimensional honeycomb lattice, has been 
synthesized^— and attracts much attention 13-15 . Sil- 
icene shares almost every remarkable property with graphene. 
Indeed, it has Dirac cones akin to graphene. It has addition- 
ally a salient feature, that is a buckled structur e) 13 ' 14 owing to 
a large ionic radius of silicon (FigJT}. Consequently, silicene 
has a relatively large spin-orbit (SO) gap of 1.55meV, which 
provides a mass to Dirac electrons. Furthermore, we may 
control experimentally the mass by applying the electric field 
E z perpendicular to the silicene sheet. As \E Z \ increases, 
the Dirac mass decreases linearly, and vanishes at the critical 
point \E Z \ = E CI , and then increases linearly 15 . 

In this paper we analyze the QH effect in silicene, which 
is the first step of exploring the intrinsic properties of silicene 
as in the case of graphene. The Hamiltonian contains the ef- 
fective spin-orbit coupling term (oc AsoX me Rashba spin- 
orbit coupling (oc Ar) and the electric field term (oc E z ) ad- 
ditionally to that of graphene. The Dirac mass depends on 
Aso and E z . There are 4-fold degenerate zero-energy states 
in graphene. Let us call the counterpart the would-be zero- 
energy states in silicene. Their degeneracy is completely re- 
solved except at the critical electric field, E z = ±E cr , due 
to the mass effect. Hall plateaux develop at v = 0, ±2 for 
E z =0,u = ±1, ±2for£ z = ±E a , and v = 0, ±1, ±2 else- 
where. Hence the critical electric field E a - can be experimen- 
tally measurable. Such a measurement will be important since 
it has been argued 15 that a transition occurs between topolog- 
ical and band insulators at this critical field without magnetic 



FIG. 1: (Color online) Illustration of the buckled honycomb lattice 
of silicene. A honeycomb lattice is distorted due to a large ionic 
radius of a silicon atom and forms a buckled structure. The A and B 
sites form two sublattices separated by a perpendicular distance 21. 
The structure generates a staggered sublattice potential in the electric 
field E z , which leads to various intriguing pheneomena. 



field. Our finding will be the first case in which we can control 
the Landau level by applying external electric field parallel to 
the magnetic field. 

This paper is composed as follows. In Section [TT] we cal- 
culate the energy levels in the silicene system under magnetic 
field as a function of E z . This allows us to determine the se- 
ries of filling factors at each E z . In Section ITlTl we elucidate 
fully the physical meaning of the energy spectrum. When we 
switch off the Rashba coupling (Ar = 0), the Hamiltonian 
is diagonalized analytically and the physical picture becomes 
manifest. Two successive Landau levels are mixed by the in- 
trinsic Zeeman effect acting between the A and B sublattaices 
of a honeycomb lattice just as in the case of graphene 1 —. The 
underlying symmetry is the supersymmetry — even in the pres- 
ence of the Dirac mass term. Each energy level is 2-fold de- 
generate in general. On the other hand, as we have mentioned, 
the four states in the would-be zero energy states are not de- 
generate: Namely, each of these states contains electrons with 
a definite spin polarization either from the K or K' point. The 
effect of the Rashba coupling (Ar ^ 0) is to remove the 2-fold 
valley degeneracy by modifying the energy level in the order 
of Ar. Furthermore, it mixes an up-spin state and a down-spin 



FIG. 2: Energy levels as a function of E z for \r 7^ 0. There exist 
two levels which yield the zero energy at E z = ±£W in the K valley, 
and also in the K' valley. (We call them the would-be zero-energy 
states.) Level crossing occurs at E z = nE\, n — 0, ±1, ±2, ■ • ■ , 
except at E z — r\E\ in the K,, valley where anticrossing takes over. 
See d3 . 1 6b with respect to E\. 



state at a certain crossing point and to turn it into an anticross- 
ing point. In conclusion, all the degeneracy is removed except 
for the points at which two levels cross. 



H. LOW-ENERGY DIRAC THEORY 

Silicene consists of a honeycomb lattice of silicon atoms 
with two sublattices made of A sites and B sites. The states 
near the Fermi energy are tt orbitals residing near the K and 
K' points at opposite corners of the hexagonal Brillouin zone. 
(We also call them the K ?) points with r\ = ±.) We take a 
silicene sheet on the xy-plane, and apply the electric field E z 
perpendicular to the plane. Due to the buckled structure the 
two sublattice planes are separated by a distance, which we 
denote by 2£ with £ = 0.23A , as illustrated in FigQ] It gener- 
ates a staggered sublattice potential oc 2£E Z between silicon 
atoms at A sites and B sites. 

We analyze the physics of electrons near the Fermi energy 
more. We employ the low-energy Dirac theory, which has 
been proved to be useful in the study of graphene^— . The 
low-energy effective Hamiltonian around the K v point reads 



FIG. 3: (Color online) Energy levels as a function of E z for Xr = 
0. All electrons are spin-up (red curves) or spin-down (blue curves) 
polarized. There exist two levels which yield the zero energy at E z = 
±E a in the K valley, and also in the K' valley. Level crossing occurs 
at E z = nE\, n — 0, ±1, ±2, • ■ • . There exists no anticrossing. 
Compare the spectrum with the one in Fig[2] 

asil 

H v = hv ¥ {k x r x - rjkyTy) + r]T z h n + £E z t Zi (2.1) 

with 

/in = -AsoCz - aX R (k y (T x - k x <7 y ) , (2.2) 

where t q is the Pauli matrix of the sublattice, = ^at = 
5.5 x 10 5 m/s is the Fermi velocity, a = 3.86A is the lattice 
constant, Aso = 3.9meV is the effective spin-orbit coupling, 
and Ar = 0.7meV is the intrinsic Rashba spin-orbit coupling. 
The two Hamiltonians H+ and H- are related through the 
time-reversal operation. See the Appendix with respect to the 
tight-binding model and the low-energy Dirac theory in sil- 
icene. 

The energy spectrum is readily derived from d2.lt as 




(2.3) 

The gap is given by 2| A^ Sa (E z ) | with 

A,., (E z ) = -?7.s,A S o + £E Z . (2.4) 
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As \E Z \ increases, the gap decreases linearly, and vanishes at 
the critical point \E Z \ = E CT with 



E„ = r]s z X S o/e = 17meV/A, 



(2.5) 



and then increases linearly. Note that the spectrum depends 
on the spin index s z and the valley index 77 only in the combi- 
nation rjs z . 

We apply a homogeneous magnetic field B — V x A = 
(0, 0, —B) with B > along the z axis to silicene. By making 
the minimal substitution, the Hamiltonian is given by 



H n = hv F (P x t x - T]P y T y ) + r)T z hu + iE z r : 



(2.6) 



with the covariant momentum Pi = ki + eAi. We introduce a 
pair of Landau-level ladder operators, 

_ e B (p x + iP y ) t _ e B (p x -iP y ) n 
V2h ' V2h ' n 



satisfying [a, a^] = 1, where l B = ^Jh/eB is the magnetic 
length. In the basis {V'At' ipBf, i>A\,i ^PsiY , the Hamiltonian 



H v reads 



' A + (E z ) hio c a 

hoj^ -A+ (E z ) 

A_ (E z ) huj c a 



\ 



• \/2ha\R a 



■ V2ha\ R ~f 
4 £ B " 



V 



i ^ XR a hu) e tf -A_ (£ z ) y 



at the K point, with u; c = \/2hvp/£B, and 
/ A_ (£7,) ftw c at -z ^ aA " at 



(2.8) 



hio c a 



V 







\ 

A_ (E z ) i V^Ar 

A + (E z ) hw c a) 

huj c a — A + (-E z ) / 



4 " 



(2.9) 

at the K' point. Here the index + of A + (E z ) implies r\s z — 
+, that is, either (s z ,r)) = (+,+) or (s z ,ri) = (-,-). 

Inspecting the Hamiltonian H v we see that the eigenstate is 
of the form 

= (< t \N) , < f +1 |JV + 1>, n^" 1 \N - 1) , \N)) 1 , 

(2.10) 

and 

* W = (< t +1 \N + l), < t I JV) , |JV) , <^ |JV - 1))* , 

(2.11) 

with 



\N) 



a tJV |0) , 



(2.12) 



which represents the state in the Ath Landau level. It is no- 
table that an energy eigenstate is a mixture of four states com- 
ing from three different Landau levels in general. The state 



(12.10b is defined for N > 1. There exists two more states 
corresponding to N = — 1 and N = in ( 12.10b . 



and 



^04. 



= (|0) ,0,0,0)% (2.15) 
*^=K t |l),t;S t |0),»^|0) J 0)*. (2.16) 

They exhaust all the eigenstates of the Hamiltonian H v . Ap- 
plying the Hamiltonian H v to these states, we can determine 
numerically the energy spectrum as a function of E z , which 
we display in Fig|2] All energy levels are nondegenerate ex- 
cept for certain isolated values of E z including E z = 0. In 
conclusion, the 4-fold degeneracy in graphene is completely 
resolved except for these isolated points, which yields the QH 
plateaux at v = 0, ±1, ±2, ±3, • • • . It is remarkable that the 
degeneracy associated with the K-K' symmetry is resolved in 
silicene. 

It is instructive to set Ar = in the Hamiltonians ( 12.81 l and 
( 12.91 l. since the Rashba coupling Ar is a small parameter com- 
pared with the others. The Hamiltonian becomes block diag- 
onal, where the spin s z is a good quantum number. Namely, 
each energy level has a definite spin polarization at the K and 
K' point separately. The resultant Hamiltonians are exactly 
solvable. We give the energy spectrum as a function of E z in 
Figf3] By comparing the spectrum for Ar 7^ and Ar = 0, it 
is found that the spectra look almost the same between them 
except for a certain value of E z , that is, E z = r\E\ in the 
valley. At this value a simple level crossing occurs for 
Ar = 0, but an anticrossing takes place for Ar 7^ 0. Though 
the spectra look quite the same, there exists actually a correc- 
tion of the energy of the order Ar in the system with Ar 7^ 0, 
which removes the degeneracy from all levels except for level 
crossing points. 

III. ENERGY SPECTRUM 

We explore the energy spectrum in detail. The silicene 
Hamiltonian dAlb is reduced to the graphene Hamiltonian if 
we set Aso = Ar = and I = 0, where there exists the 4-fold 
degenerate zero-energy states. Let us call them the would-be 
zero modes in silicene. We wish to study how the degeneracy 
is removed. We discuss the states ( 12.13b . ( 12.14b . and ( 12.10b in 
this order. 

First of all, we note that the leading order of the Hamilto- 
nian H v takes a simple form in the limit E z — > ±00, 

H v = £E Z [diag.(l, -1, 1, -1) + 0(1/E Z )\. (3.1) 

Consequently, there are four eigenstates in this limit such that 

(3.2a) 

(3.2b) 

(3.2c) 

(3.2d) 



(0,|0),0,0) > (2-13) 
« T |0),^ T |1),0,4J0))*, (2.14) 



(1) 


u A-\ — 


1, 


a Bt 




1 




= 0, 


(ii) 


u N+1 


= 1, 






1 




= 0, 


(iii) 




= 1, 


u A-\ 


— a B ^ 


1 




= 0, 


(iv) 


u N Bi = 


1, 


v N - 


u N+1 






= 0. 
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FIG. 4: (Color online) Energy levels as a function of E z for the state 
$2^ and There are four levels (solid) from the K point and four 
levels (dotted) for the K' point. The red lines are the would-be zero 
mode made entirely of up spins. On the other hand, the blue thin lines 
are the would-be zero mode made of entirely of down spin if Xr = 0. 
They anticross with other levels with up spin at E z — ±E\, and 
become mixed states when \r 7^ 0. It is remarkable that there exists 
no degeneracy for < \E Z \ < E a . Namely, both the spin and valley 
symmeties are explicitly broken, where the QH plateaux appear at 
v — 0, ±1, ±2. The symbol such as NA^ in the left (right) column 
indicates that the state is dominated by up-spin electrons at the A site 
coming from the iVfh Landau level as E z — > — 00 (E z — > +00). 



It turns out that we can label each energy level by the two 
asymptotic states in the limit E z — > ±00. For instance, the 
label [TV Af, NBl] implies that we have = 1 as E z -> 
= 1 as E z — > +00. 



-00 and 1 



A. Would-be Zero-Energy States 



For the state ( |2.13t it is trivial to solve the eigenvalue prob- 
lem, 



H s z _ , A,, a , (E z ) fuv c a 

hu c a^ -A VSz (E z ) 

H s z _ 1 A riSz (Eg) hu c a^ 

fku c a -A VSz (E z ) 



(3.6a) 
(3.6b) 



It is easy to diagonalize H** analytically, where the spin s z is 
a good quantum number. 

Now it is trivial to obtain one of the eigenstates ( 12.141 ) at the 
K point such that u° A ^ — u l B ^ = u A ^ = and u B ^ = 1. The 
state is labelled by [OB^OBJ,], as indicated by a blue thin line 
in Fig|4] It leads to the emergence of the zero-energy state 
at the critical point E z = —E a . When we include the effect 
of Ar 7^ 0, however, an anticrossing takes place between the 
level and another level at E z = E\ as indicated by a blue 
line, and the resulting state becomes [OB^IBJ,]. We discuss 
the anticrossing later. 

Similarly we can treat H_. The level is indicated by a 
blue dotted thin line in Fig0] It leads to the emergence of the 
zero-energy state at the critical point E z = E a -. As an effect 
of Ar 0, an anticrossing takes place with another level at 
E z = -E x . 

We have discussed the four states coming from the K and 
K' points. They are the would-be zero-energy modes in sil- 
icene, which become degenerate zero-energy states in the 
limit Aso = Ar = 0. The degeneracy is completely resolved 
except at the critical electric field, E z = ±E CI , due to the 
mass effect. Hall plateaux develop at v = 0, ±2 for E z = 0, 
v = ±1, ±2 for Eg = ±E a , and v = 0, ±1, ±2 elsewhere. 
Hence the critical electric field E a - can be determined experi- 
mentally. 

There exists two more states of the type ( 12.141 ), which 
are indicated by black curves with the label [lBt,0Af] and 
[0At,0BJ,] in Fig|U We discuss them in the succeeding sec- 
tion. 



H W 01 " = e T * 01 ' 



where the eigenvalue is exactly given by 



4 = - V A V (Eg) , 



B. Nonzero Modes 



(3.3) 



(3.4) 



with ( 12.41 ). This is a reminiscence of the gap energy in the sys- 
tem without magnetic field. In particular it implies the emer- 
gence of the zero-energy state (ej. = 0) at the critical point 
Eg = r\E a . We illustrate the energy levels by red lines (solid 
line from the K point and dotted line from the K' point) in 
Fig|4] The label of this state is [0Bt,0Bt] at the K point and 
[0At,0At] at the K' point according to our convention. 

There are three states of the type ( 12.141 ). It is not easy to 
solve them analytically. However, when we set Ar = 0, the 
zero-th order Hamiltonian becomes block diagonal, 



H (0) = (^ 
\ Hjj 



(3.5) 



We go on to analyze the eigenstate ( 12.101 ), i.e., 

= (< t \N) , u^ 1 \N + 1), u^ 1 \N 1) 



(3.7) 

There are four eigenstates for each N. We set Ar = 0, where 
the Hamiltonian becomes block diagonal as in ( 13.5b . The 
Hamiltonian ( 13.61 ) is a Dirac Hamiltonian with spin dependent 
mass. The four eigenstates read as follows. 

With respect to the up-spin sector, the eigenvalue problems 
are solved as 



Hi 



where 



u^\N + 



1) 



II 



1) 



it 



N 



=[A + (Eg)±el}/C, 



=hoj c y/N 

N-l N 
Al — u Bi 



1/C, 



II 



0. 



, (3.8) 

(3.9a) 
(3.9b) 
(3.9c) 
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^ 1B4- 



FIG. 5: Energy levels as a function of E z for Xr ^ comming from 
the state ty]).. Each level is dominated by up or down-spin electrons 
at A or B sites in the limit E z — > ±oo, as indicated, but they mix 
near the anticrossing point (Ez — E\). 



with C the normalization constant. The eigenvalue is 



t 



with 



e N = yJ(?uj c ) 2 N + AZ 3z (E z ). 



(3.10) 



(3.11) 



They are labelled by [iVAt,(iV+l)Bt] and [(iV+l)Bt,iVAt]. 
See Fig[5]for the case of N = 1. 

With respect to the down-spin sector, the eigenvalue prob- 
lems are solved as 



u^\N-l) 
u^\N) 



±e 



N 



u^\N-l) 
u^\N) 



where 



=hjj c \J N /C, 
with C the normalization constant. The eigenvalue is 



JV _N+1 
a A\ ~B-\ 

„,N 



0. 



=[-A_ (E z )±e l N ]/C, 



J. 

-N 



(3.12) 

(3.13a) 
(3.13b) 
(3.13c) 

(3.14) 



with dXTTl They are labelled by [(N - l)A|,iVB4.] and 
[NBUN ~ 1)A|]. See Fig|5]for the case of N = 1. We 
can similarly discuss the K' point. 

In the previous subsection we have studied ^f^_, where 
we have pointed out that there are two levels labelled by 
[lBt,0Af] and [0At,0B!], which are nonzero modes (see 
FigHJl. The are actually the states just studied in ( 13.9b by 
choosing N = 0. 

We illustrate the energy levels in Fig[3] One level crosses 
with other levels. The crossing point is given by en — £N' in 
general, which is solved as 



E, 



(N-N')Ex, 



(3.15) 



E x = {hw c f /2X S0 £. 



(3.16) 



The spectrum with Ar = is found to be almost identical to 
that of the full theory (Ar ^ 0) in Fig|2]except for the vicinity 
of E z — r\E\ in the K,, valley, where an anticrossing takes 
place. 



C. Landau-Level Mixing and Anticrossing 

We make a further study of the physical meaning of the en- 
ergy levels. It is a prominent feature that each energy eigen- 
state is a coherent superposition of states belonging to differ- 
ent Landau levels. We explain the mechanism how such a 
Landau-level mixing occurs by setting A R = and then dis- 
cuss the anticrossing by taking the effect of Ar 7^ 0. 

When Ar = 0, the Hamiltonian is block diagonal as in 
(13.5b . where the block elements are 



with 



a 

at 



Q- 



a f 
a 



(3.17) 



(3.18) 



To reveal the intrinsic structure of the energy spectrum, we 
note that 



(3.19) 



and explore the Hamiltonian Hp = Q±Q±. We explicitly 
discuss 



h; 



oat 
a)a 



but a similar result follows also for H p . We set 





at 



it _ 



a 




(3.20) 



(3.21) 



It is easy to see 

H+ = {Q,Q^}, 



[H+,Q]=0. (3.22) 



It defines a SUSY algebra, with Q the supercharge . It in- 
terchanges a fermion state and a boson state within a single 
energy multiplet called a supermultiplet. The SUSY is the 
underlying symmetry of the QH effect in silicene just as in 
graphensii 

To construct the fermion and boson sectors explicitly we 
analyze the Hamiltonian ( 13.201 ) more in detail by rewriting it 

as 



where 





1 



a/ a + crc, 



at - 



1 




(3.23) 



(3.24) 
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(a) Energy levels 



Fermion 
Sector (A) 



Boson 
Sector (B) 



(b) Landau and Energy levels 

Landau Energy Landau 
level (A) level level (B) 




N=3 



7" 

" a 2 ft$ ^ 



N=3 



ft 71 4 
ft /I $ 



FIG. 6: (a) The illustaration of the SUSY structure. The energy spec- 
tra of the fermion sector (A site) and the boson sector (B site) are 
related by the supercharge Q. The two states on the same horizon- 
tal line have the same energy, making a supermultiplet, except for 
the ground state (so) in the boson sector, (b) The mechanism of the 
Landau-level mixing. It is illustrated how a Dirac electron with the 
energy e n comes from the ATth Landau level of the A site and also 
from the (N + l)th Landau level of the B site in the Dirac K valley. 
In this figure electrons are filled up to the level £2- 



Since {c, c} = 0, and {c, c 1 } = 1, c and are the fermion an- 
nihilation and creation operators with the ground state \0}b = 
(0, 1)* and the one-fermion state 10}^ = (1,0)*, 

c|0} B = 0, tf\0} B = \0} A , c|0} a = |0} b . (3.25) 

The boson sector \N}b and the fermion sector \N}a are con- 
structed by operating the boson creation operator a 1 ' to the 
ground states in each sectors, 



\N} A =- 
\N} B =- 



1 



fNl 
1 



it* 



\0}a = 
\0}b = 



\N) 




\N) 



The energy of the states is 

H+\N} A =(N + 1)\N} A , 

H+\N} B =N\N} b . 
The two sectors are interchanged by the supercharge, 

Q\N} A =&*c\N}a = VWTT\N + 1} B , 
Q^\N + 1} B =a$\N + 1} B = VW+1\N} A . 



(3.26) 



(3.27) 



(3.28) 
(3.29) 

(3.30) 
(3.31) 



The two states \N}a and \N + 1}b make a supermultiplet. 
We may identify \N}a as the iVth Landau level at the A site, 
and \N + 1} B as the (iV + l)th Landau level at the B site. 
They are degenerate, as illustrated in Fig Ob). An eigenstate 
of the Hamiltonian fff is a coherent superposition of them as 

in ( I3.9l l. Similarly we can analyze the Hamiltonian H\. 

The physical reason of the degeneracy is understood as 
follow a 16 ' 17 . We rewrite ( 13 . 1 9b in the following form, 



{-ihV + eAf + V ehBT 3 + A* Sz (E. 



This is essentially the Pauli Hamiltonian. The first term is the 
kinetic term, while the second term is the Zeeman term cre- 
ating the energy gap between the A site and B site. The Lan- 
dau level is generated by electrons making cyclotron motion, 
which yields the energy (N + ^)(hu! c ) 2 in the TVth Landau 
level. On the other hand, there exists the intrinsic Zeeman en- 
ergy ±^(fno c ) 2 between the A and B sublattices 16 , and their 
sum is either N(huj c ) 2 or (N + l)(ftw c ) 2 . Consequently, the 
energy of massive Dirac electrons is either ie^v or ±ejv+i, as 
given by (13. lit . Namely, Dirac electrons with the energy ejv 
come from the AAh Landau level of the A site and also from 
the (N + l)th Landau level of the B site in the Dirac K valley, 
as illustrated in Fig|6jb). 

Finally we discuss the effect of Ar 7^ 0. The Hamilto- 
nian is no longer block diagonal, where off-diagonal block 
elements mix up-spin states of Hi. and down-spin states of 
Hl_. The mixing occurs solely within the four states de- 
scribed by J2. 10b . It is enough to consider a mixing between 
positive-energy states labelled by [(N + l)Bt,(N-l)A|] and 
[NBl,(N - l)Al]: See Fig|5]for the instance of = 1. They 
have the energy e^ N — En+i and — £jv, respectively, as 
we have shown in d3.9t and (13.13b . When Ar = 0, they cross 
at Ejf — ej^, which implies En+i = £n with (13.11b . Solving 
this equation we find 



E z = r]Ex 



(3.33) 



in the K v valley, with E\ being given by (13. 16b . When Ar 7^ 
0, the Rashba interaction operates between these two levels. 
This is a typical quantum mechanical problem. As a result 
of the interaction, an anticrossing takes place with the energy 
gap 



2V2ha , 
oe = — Ar 



(3.34) 



(3.32) 



at the point ( 13.331 ). There are many level crossings, but anti- 
crossings occur only at this point. 



IV. CONCLUSIONS 

We have studied the QH effect in silicene, which is a mono- 
layer of silicon atoms replacing carbon atoms in graphene. 
The system is described by a massive Dirac theory due to 
a large spin-orbit couplings. We have calculated the energy 
spectrum under magnetic field, which exhibits an intriguing 
feature as a function of external electric field as in Fig|2] We 
have unveiled a rich physical picture behind it. Our finding 
will be the first system in which the Landau level is controlled 
by applying external electric field parallel to the magnetic 
field. 

I am very much grateful to N. Nagaosa for many fruitful 
discussions on the subject. This work was supported in part 
by Grants-in-Aid for Scientific Research from the Ministry of 
Education, Science, Sports and Culture No. 22740196. 
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Appendix A: Tight-Binding Model 

In this appendix we derive the low-energy Dirac Hamilto- 
nian ( 12.lt in silicene. The silicene system is described by the 
four-band second-nearest-neighbor tight binding modeUi, 



first and the second terms constitute the Kane-Mele model 23 
proposed to demonstrate the QSH effect in graphene. 

Taking the continuum limit is well known except for the 
Rashba term. With respect to the Rashba term, we take the six 



vectors <£?„■ as 



«*,J>>«/3 



((hi))'** 



(Al) 



The first term represents the usual nearest-neighbor hopping 
on the honeycomb lattice with the transfer energy t — 1.6eV, 
where the sum is taken over all pairs of the nearest- 

neighboring sites, and the operator cj a creates an electron 
with spin polarization a at site i. The second term rep- 
resents the effective SO coupling with Aso = 3.9meV, 
where a — (a x ,a v ,<j z ) is the Pauli matrix of spin, i/y = 

di x dj with di and dj the two nearest bonds 



(di xdj^j I 



connecting the next-nearest neighbors, and the sum is taken 
over all pairs of the second-nearest-neighboring sites. 

The third term represents the intrinsic Rashba SO coupling 
with Ar = 0.7meV, where fUj — ±1 for the A (B) site, and 



dijj 



tential term, where r/i — ±1 for the A (B) site. Note that the 



The forth term is the staggered sublattice po- 



4 = 



/ nn . 7rn\ 

[cos— ,sm— J, (A2) 
with n = 0, 1, • • • ,5. Then we obtain 

E/. it)\ z r, • ak x . V3ak y 

I a x d-j I = 2ma x cos — — sin — — - 



2ai 



^-<Ty \ sin ak x + sin — — cos 



ak x ^/?>ak v 
— cos — — — 
2 2 



(A3) 



We have taken the K and K' points as 



*=l£A «•-(-£,» I. .A4, 



We may expand it as 



. z iai 



^2 (ff x ^ = ^ (a y k x - a x k y ) . (A5) 



This is the low-energy continuum version of the Rashba term 
in the Dirac Hamiltonian (12. U . 



1 A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. Novoselov and 
A.K. Geim , Rev. Mod. Phys. 81 (2009) 109. 

2 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, 
S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov: Science 306 

(2004) 666. 

3 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Kat- 
snelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov: Nature 
438 (2005) 197. 

4 Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim: Nature 438 

(2005) 201. 

5 K. S. Novoselov, E. McCann, S. V. Morozov, V. I. Fal'ko, M. I. 
Katsnelson, U. Zeitler, D. Jiang, F. Schendini, and A. K. Geim: 
Nat. Phys. 2 (2006) 177. 

6 J. Alicea and M. P. A. Fisher: Phys. Rev. B 74 (2006) 075422. 

7 K. Nomura and A. H. MacDonald: Phys. Rev. Lett. 96 (2006) 
256602. 

8 M. Goerbig, R. Moessner, and B. Doucot: Phys. Rev. B 74 (2006) 
161407. 

9 M.O. Goerbig, Rev. Mod. Phys. 83 (2011) 1193, and references 
therein. 

10 B. Lalmi, H. Oughaddou, H. Enriquez, A. Kara, S. Vizzini, B. 
Ealet, and B. Aufray, Appl. Phys. Lett. 97 (2010) 223109. 

11 P.E. Padova, C. Quaresima, C. Ottaviani, P.M. Sheverdyaeva, P. 



Moras, C. Carbone, D. Topwal, B. Olivieri, A. Kara, H. Oughad- 
dou, B. Aufray, and G.L. Lay, Appl. Phys. Lett. 96 (2010) 261905. 

12 B. Aufray A. Vizzini, H. Oughaddou, C. Lndri, B. Ealet, and G.L. 
Lay, Appl. Phys. Lett. 96 (2010) 183102. 

13 Gian G. Guzman- Verri and L. C. Lew Yan Voon, Phys. Rev. B 76 
(2007)075131. 

14 C.-C. Liu, W. Feng, and Y. Yao, Phys. Rev. Lett. 107 (2011) 
076802; C.-C. Liu, H. Jiang, and Y. Yao, Phys. Rev. B, 84 (2011) 
195430. 

15 M. Ezawa, cond-mat/arXiv:1201.3687. 

16 M. Ezawa, J. Phys. Soc. Jpn 76 (2007) 094701. 

17 M. Ezawa, Phys. Let. A 372 (2008) 924. 

18 J.C. Slonczewski and PR. Weiss, Phys. Rev. 109 (1958) 272. 

19 G.W. Semenoff, Phys. Rev. Lett. 53 (1984) 2449. 

20 H. Ajiki and T. Ando, J. Phys. Soc. Jpn. 62 (1993) 1255. 

21 T. Ando, Y. Zheng, and H. Suzuura, Microelectron. Eng. 63 
(2002) 167. 

22 E. Witten, Nucl. Phys. B 185 (1981) 513; E. Witten, Nucl. Phys. 
B 202 (1982) 253. 

23 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005); 
ibid 95, 146802 (2005). 



